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SUMMARY 


For initial stagnation pressures up to 1,000 atmospheres and stagna- 
tion enthalpies up to 10,000 Btu per pound, nozzle-flow properties for 
equilibrium air have been computed and plotted on charts . The work of 
NASA TN D -693 has "been extended to include flow properties for closer 
intervals of specified stagnation enthalpies . Properties which have been 
charted as a function of Mach number are as follows: temperature, 

pressure, density, velocity, area ratio, dynamic pressure, Reynolds number 
isentropic exponent, and molecular weight ratio. Ratios of temperature, 
pressure, and density across normal shock waves are also charted, and • 
weight-flow rate is plotted as a function of stagnation enthalpy. 


INTRODUCTION 


In order to design hypervelocity nozzles and to define the test 
conditions, various air-flow properties must be determined. Several 
theoretical investigations of air in equilibrium, frozen, and nonequilib- 
rium flows have been reported (e.g., refs. 1 to 4) . In a recent study 
(ref. 5 ) charts for equilibrium and frozen air flows were computed for 
specified stagnation pressures of 6 . 8 , 68 , and 680 atmospheres and stagna- 
tion enthalpies of 125, 2,000, 4,500, 8,000, 12,000 and 24,500 Btu per 
pound. In the present study the work of reference 5 has been extended to 
include equilibrium air-flow properties for stagnation pressures up to 

1.000 atmospheres and for many additional stagnation enthalpies below 

10.000 Btu per pound. Also, in contrast to reference 5, the flow proper- 
ties have been plotted as a function of Mach number instead of area ratio 
for convenience of use in wind-tunnel operation. 

The purpose of this report is to present for specified stagnation 
pressures and enthalpies the following properties as a function of Mach 
number: temperature, pressure, density, velocity, area ratio, dynamic 

pressure, Reynolds number, isentropic exponent, and molecular weight ratip 
Also included are temperature, pressure, and density ratios across normal 
shock waves . Weight-flow rate is plotted as a function of stagnation 
enthalpy. Because in an actual nozzle the flow may not be in complete 
equilibrium, care should be exercised in the application of these charts . 
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NOTATIONS AND CONSTANTS 

nozzle cross-section area, ft 2 
isentropic speed of sound, ft/ sec 

enthalpy (with datum of h = 0 for molecular gas at 0° K] 
enthalpy, dimensionless 
characteristic length, ft 
Mach number 

molecular weight of mixture, lh/mole 

molecular weight of undissociated air, 28.966 lh/mole 
pressure, atmospheres unless specified 
dynamic pressure (l/2 pu 2 ) , lb /ft 2 
Reynolds number , upL/ p 

gas constant, 1.987 Btu/mole-°R or 0.7302 ft 3 -atm/mole- . 
33.73 Btu/lh for air 
entropy, Btu/( initial mole) (°R) 
entropy, dimensionless 
temperature, °R or °K as specified 
reference temperature, ^92° R = 273° K 
velocity, ft/sec 
weight-flow rate, lh/sec 
molecular weight ratio, mo/m 
isentropic exponent, 
density, slug/ft 3 

reference density, 0.00251 slug/ft 3 for air 


Ta v \ 

In P/s 


, Btu/lh 


v 
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p coefficient of viscosity, lb-sec/ft 2 

( )* sonic point 

Subscripts 

0 reference condition 

t reservoir condition 

1 condition in front of normal shock wave 

2 condition behind normal shock wave 


THERMODYNAMIC PROPERTIES 


For ready reference and convenience, thermodynamic properties of air 
for temperatures to 10,000° K and pressures to 10 3 atmospheres are plotted 
in figure 1. These plots give the variation of molecular weight ratio Z 
(fig. 1(a)), enthalpy h/RT 0 (fig. 1(b)), and entropy S/R (fig. 1(c)) with 
temperature for specified pressures. The reader is referred to the previ- 
ous section of the report for values of the constants RT 0 and R. For 
temperatures to 3,000° K and pressures from 10 -2 to 10 2 atmospheres, the 
properties were obtained from reference 6. For temperatures from 2,000° K 
to 10,000° K and pressures from 10" 5 to 103 atmospheres, the properties 
were obtained from reference 7* The properties from references 6 and 7 
were in excellent agreement in the overlapping range of temperatures and 
pressures (2,000° K to 3,000° K and 10“ 2 to 10 2 atmospheres). In addition 
to information from these sources, some of the curves at low pressures and 
temperatures (p as low as 10~ 6 atmosphere and T as low as 50° K) were 
determined based on the perfect gas law and are identified by dashed lines. 
For more extensive air charts for temperatures to 15,000° K, the reader is 
referred to references 8 and 9; and for air properties at temperatures to 
100,000° K, the reader is referred to reference 10. 

Vapor pressures as a function of temperature for air, oxygen, and 
nitrogen are presented in figure 2. To estimate the vapor pressures for 
air a perfect liquid or solid solution of oxygen and nitrogen was assumed 
in equilibrium with a 20-percent oxygen, 80-percent nitrogen vapor. The 
vapor-pressure data for oxygen and nitrogen were obtained from references 
6, 11, and 12. The approximate saturated vapor line for air is also 
indicated on the entropy graph in figure l(c) . 



CALCULATION OF NOZZLE FLOW PROPERTIES 
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The nozzle-flow properties as a function of Mach number were desired 
for use in wind-tunnel testing. For convenience in the computational 
procedure , however, the flow properties were computed as a function of the 
ratio of static to stagnation pressure, with the corresponding Mach number 
being computed. Stagnation conditions p t and ht at the reservoir were 
specified for each case considered. In order to compute the area ratio, 
Mach n umb er, and Reynolds number, the density was determined from the 
thermal equation of state 


p = pmo/ZRT 


( 1 ) 


For p in slugs per cubic foot , p 


in atmospheres, and T in °R, 


p = 1.23^ p/ZT 


( 2 ) 


where 

P - P t (p/Pt) 


An isentropic expansion of the gas from the reservoir to downstream 
stations in the nozzle was assumed, and the molecular weight ratio Z, 
the temperature T, and the enthalpy h = (h/RT 0 )RT 0 were determined from 
the ther mal plots (fig. l) . The velocity was then computed from the energy 
relation 


u = 223.6 - h 


(3) 


The area ratio was calculated from the continuity equation 


A/A* = p*u*/pu (4) 

where the mass-flow rate at the throat (w/A* = p*u*) was determined as the 
maximum pu for given stagnation conditions p^ and h^. The Mach number 
was computed from 

M = § *> U ■■■ ■ (5) 

Vtzrt 


For T in °R the speed of sound a in feet per second was given by 


VJ1 f— 1 o\ > 
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a = 41.43 n/tZT 


( 6 ) 


where the isentropic exponent y was determined from 


7 = 


The Reynolds number was calculated from 

( 8 ) 

LPt W?t 

with the coefficient of viscosity p computed hy the method of refer- 
ence 13. In reference 13 the effects of molecular dissociation and ioni- 
zation on p were considered. 

An iterative procedure was used to obtain equilibrium temperature 
ratios , pressure ratios, and density ratios across a normal shock wave. 
This procedure is outlined in references 3 and Taut it is repeated here 
for clarity and convenience. Since the Mach number Mi in front of the 
shock was known, the density ratio, p / p i , was assumed, and the pressure 
behind the shock was calculated from the expression 

P 2 = P 1 + ( p l U l) (u x ) (l - p^) (9) 

which resulted from combining the momentum and continuity equations . For 
the pressure in atmospheres 

P 2 = p i + 2^7 (( W (u i ) ( P ' §) (10) 

The local enthalpy behind the shock, tu, was determined by the energy 
equation 



h 2 = h x + ( l/2)u % 2 - (l/2)u2 2 


(ID 


For h 2 in Btu per pound 


h 2 = hi + 


50,060 


U1 2 



( 12 ) 


With p 2 and h 2 known, T 2 and Z 2 were obtained from the thermal plots 



6 


(fig. l) . The density "behind the normal shock, p 2 , was then computed from 
equation (l), and p 2 /p was compared with the assumed value. This process 
was repeated until the assumed and computed values of p 2 / p 1 were in 
agreement. The stagnation pressure behind the normal shock was then deter- 
mined from the thermal plots with an isentropic compression assumed from 
h 2 and p 2 to h t - h t . 


RESULTS 


The procedure outlined in the previous section was used to compute 
nozzle properties for equilibrium flow for various initial stagnation 
pressures and enthalpies . All calculations were made on an IBM 7090 
system, the thermodynamic properties of air (from the references previously 
cited) hav ing been put on magnetic tape . The nozzle-flow properties 
determined were temperature, pressure, density, velocity, area ratio, 

Mach number, dynamic pressure, Reynolds number, isentropic exponent, 
molecular weight ratio, and weight flow. Ratios of temperature , pressure, 
and density across normal shock waves were also determined. The initial 
stagnation pressures for the calculations were 1, 10, 100, and 1,000 
atmospheres, and the initial stagnation enthalpies included values from 
4-00 to 10,000 Btu per pound. 

Except for weight-flow rate which was plotted as a function of stagna- 
tion enthalpy, al 1 properties were plotted as a function of Mach number 
and are presented on charts which are indexed on page 17. The reader 
should note that, except on the weight-flow chart (chart l4) , symbols are 
used on the charts only to identify curves, not to indicate the many points 
computed for each curve. Perfect-air properties (from ref. l4) are 
presented on most of the charts for comparison with the real-air 
properties . 


Ames Research Center 

Rational Aeronautics and Space Administration 
Moffett Field, Calif., June 27, 19^2 
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Figure 1.- Thermodynamic properties of air. 
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Figure 1.- Concluded. 
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Figure 2.- Vapor pressures of air, oxygen, and nitrogen. 
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Chart 2.- Variation of temperature ratio across a normal shock wave with 

Mach number. 
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Chart 3.- Variation of the ratio of static to total pressure with Mach 

number . 



Chart 3*- Continued. 


\J1 i— 1 ON > 



lO ' 5 


I0' 6 


lO ’ 7 


8 9 10 II 12 13 14 15 16 17 18 19 20 21 22 

M 

(b) = 10 atm - Concluded. 


ro 

vo 


vji H av > 


Chart 3*- Continued 


r 



(c) p t = 100 atm 


Chart 3 Continued. 


J 



Chart 3*- Continued. 


uu 

H 



icr 


10 " 


p 

Pt 


10 “ 


IO- ( 


(d) p t = 1000 atm 


U) 

ro 


Chart 3*- Continued. 


VJl H Q\> 


u \ ^ 


p 

p 4 



10' 


10 " 


p 

Pt 


10 " 


10 " 


do 

> 


(d) = 1000 atm - Concluded. 


Chart 3.- Concluded. 


uj 

uo 


/ — 



Chart k.- Variation of pressure ratio across a normal shock wave with 

Mach number. 
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Chart 5.- Variation of total pressure ratio across a normal shock wave 

with Mach number. 
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Chart 6.- Variation of density parameter with Mach number; 
p Q = 0.00251 slug/ft 3 . 
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Chart 7*- Variation of density ratio across a normal shock with Mach 

number. 
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Chart 8.- Continued. 
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Chart 9 *- Variation of area ratio with Mach number 
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Chart 10.- Variation of dynamic pressure with Mach number. 
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(b) p t = 10 atm 


Chart 10.- Continued 
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(c) p f =100 atm 
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Chart 10.- Continued. 
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(d) p t =1000 atm 
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Chart 10.- Concluded 
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Chart 11.- Variation of Reynolds number with Mach number,. 
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Chart 11 
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Chart 11.- Continued. 
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Chart 11.- Continued 
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Chart 11.- Continued 
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(d) p t =1000 atm 




Chart 11.- Concluded 
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(d) p f =1000 atm 

Chart 12.- Variation of isentropic exponent with Mach number 
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Chart 13*- Variation, of molecular weight ratio with Mach number. 
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Chart 1^.- Variation of weight-flow rate at throat with stagnation enthalpy. 
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